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Abstract: 

The minimal conductivity of graphene is a quantity measured in the DC limit. It is shown, using the 
Kubo formula, that the actual value of the minimal conductivity is sensitive to the order in which certain 
limits are taken. If the DC limit is taken before the integration over energies is performed, the minimal 
conductivity of graphene is 4/tt (in units of e 2 /h) and it is tt/2 in the reverse order. The value tt is 
obtained if weak disorder is included via a small frequency-dependent selfenergy. In the high-frequency 
limit the minimal conductivity approaches tt/2 and drops to zero if the frequency exceeds the cut-off 
energy of the particles. 



Introduction. The conductivity <r w of graphene varies with the density of quasiparticles almost 
linearly with a minimal value cr mm « 4e 2 /h p] [2]. In terms of theoretical calculations, there has been 
some confusion about the actual value of the minimal conductivity cr mm . This confusion is twofold: 
one originates from the experimentally observed value that is roughly three times as big as most of the 
calculated values, the other one is related to the theoretical calculations that produced different values 
of a mm (calculated per spin and per valley): 

°f u = lj lasiEiEimEiiiiiini nn, ^ min = ^ umm, ^ in = ^ m- a) 

erf 111 was obtained from the Kubo formula [3[ [H [6l [7l [10l [11] as well as from the Landauer formula 
[3 [HI Qj] , whereas a™™ were obtained from the Kubo formula only. All these results were calculated 
near the ballistic regime of the quasiparticles. The possibility to reach the experimentally observed 
values of the minimal conductivity by including long-range disorder due to charged impurities was also 
discussed recently [2]. The latter will not be considered in the subsequent discussion. Instead, it shall 
be explained that all the results in Eq. fl} can be obtained from the standard Kubo formula of nearly 
ballistic quasiparticles by taking limits in different order. When a non-zero temperature T is considered, 
the conductivity is a function cr mm (w/T, rj/T), for frequency u> and scattering rate rj. 

Frequency-dependent conductivity. The conductivity a^ v is given by the Kubo formula as a response 
to an external field with frequency to. Here the representation given in Ref. [13) is used 

^,=4/ fTr{{H,r,]5(H-e'){H,r„}S(H-e)\ M£1_JM dede' , (2) 

n J J I J e - e + uj — la e — e 

where //3(e) = 1/(1 + exp(/3e)) is the Fermi function at temperature T = l/(/cs/3)- For the minimal 
conductivity only the real part of the diagonal conductivity a' uu = Re(a vv ) is of interest. After taking 
the limit a — * 0, the e' integration can be performed and gives 

o-l - 4 / H t g > r ^ H - e - ^ r ^ H 4 — + "i " fp{£) de - w 

In the zero-temperature limit (3 — > 00 this becomes 

o' vv = -^-\' Tr{[H,r„]5(H - u/2 - e)[H,r v ]5{H + tu/2-e)\de. (4) 

11 w J-u/2 L ' 
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The minimal conductivity cr mm is obtained by taking the limit uj — > of a' vv . It is tempting to ignore 
the uj dependence of the integrand and replace the right-hand side by the integrand at e — uj — 0. It will 
be shown subsequently that this does not agree with the result, when we perform the energy integration 
first and take the limit uj — > later. 

Dirac fermions. The Hamiltonian of Dirac fermions in 2D with wavevector (fci, fo) 

H = o\ki + a 2 k 2 (5) 

describes the low-energy quasiparticles in graphene. Oj (j = 1, 2, 3) are Pauli matrices. H can be diago- 
nalized as diag(k 1 —k) with k — \/k\ + k 2 . The current operator transforms under Fourier transformation 

as 

Jfl = -te[H,r„\ ~^ e g^- 

This means that the current operators for the Hamiltonian H in Eq. ([5]) is a 2 x 2 matrix with vanishing 
diagonal elements. The representation of j% in terms of energy eigenstates reads 

e ( k 2 iki 



h= k[-ih -fcj- (6) 

Thus, the current j 2 does not depend on k but only on the polar angle. The trace term in the conductivity 
of Eq. (gj) reads 

T(e) = -Tr{[H, r^]5(H - uj/2 - e)[H, r^]8{H + uj/2 - e)} 

^{^(H-^-^H + ^-e)}^, (7) 

where Tr 2 is the trace with respect to 2 x 2 matrices. After diagonalizing H this becomes together with 
the current in Eq. © 

f k 2 rlkitlkr, 

T(e) = / -i[<5(fc + uj/2 + e)5(k + uj/2 - e) + S(k - uj/2 - e)S(k - uj/2 + e)]^f§- 
J k (^ 7r ) 

k 2 rlkiflkn 

M5{k - uj/2 - e)S(k + oj/2 - e) + 5{k + uj/2 + e)S(k - uj/2 + e)]-^-# (8) 

k (27Tj 

which is a symmetric function with respect to e. 

Now a soft Dirac delta function S v (x) is considered with 



v[X> ~ nx 2 +r, 2 2nr 
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x + irj x — ir\ 



(9) 



The parameter r\ (a scattering rate) can be understood as the imaginary part of the selfenergy, created, 
for instance, by random fluctuations due to disorder [4] . With the energy cut-off A for the Dirac fermions 
and rj ~ 0, the integral of the double product of soft Dirac delta functions reads 

S v {k ~ a)S v (k - b)kdk - (a + b)5 v (a ~ 6)i [9(A - a) + 9(a) + 9(A - b) + 9(6) - 2] 
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11 1 [6(A - 6) + 6(6) - 9(A - a) - 6(a)] . 



a — 6 2tt 

Returning to Eq. ([8]), we restrict the variable e to —uj/2 < e < w/2, since for low temperature the term 

sinh(Buj/2) 

M e + uj/2) - M e - uj/2) = - ^^/^ 
in the conductivity is exponentially small for |e| > oj/2. If it is further assumed that 77,0; A we obtain 

T W ~ TTv (7 W + — ) Q ( A - w / 2 )> ( 10 ) 

\2tx) 1 V 4 TTUJ/ 
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where the prefactor tt is a result of the angular integration of k 2 jk 2 . The first term describes interband 
scattering (i.e. scattering between states with different energies ±e) and the second term intraband 
scattering (i.e. scattering between states with the same energy e or — e). The intraband scattering term 
increases linearly with the scattering rate ?y, in contrast to the 77-independent interband scattering. The 
frequency dependence is also different for the two types of scattering: the interband term increases with 
uj, whereas the intraband term decreases. 

The temperature-dependent conductivity can be calculated from Eqs. ([3]) and (J?) as 



hi oj 



, e 2 77 r /2 sinh(/W2) , , , 



Thus Eq. (JTUJ) implies 



^ 8h lJpy 1 ' Jpy ' n 2huj 2 J_ u/2 cosh(/3e) + cosh(/?w/2) 
for uj < 2A and a vanishing conductivity for uj > 2A. The integral in the second term gives 
f r /2 sinh(/3w/2) 4 



f r jA sinh (0w/2) , 4 , r , %rn 

- / y 1 n \ , 77 ,^ de = -^arctanh tanh 2 (/3w/4) 

w i_ w /2 cosh(/3e) +cosh /3w/2) ;J 



Moreover, the relation 



can be used to get 



arctanh(x) = — log 



; Tre 2 , , - , ,, e 2 3r> , 
4^-tanh( /3 ./4) + - I |l I log 



1 + tanh 2 (/3w/4) 
1 -tanh 2 (/3w/4) 



(13) 



This is the main result of this paper, ft shows that the conductivity depends on two parameters, (3uj 
and /3rj. Experimentally interesting is the case where (3rj is fixed and [3uj is varied (cf. Fig. 1). This is 
motivated by two facts. The first one is related to the origin of ~q (the scattering rate or inverse scattering 
time). It grows with increasing disorder. An important source of disorder in graphene are "ripples" in 
the carbon sheet [TUBS] which are created by thermal fluctuations. Therefore, a simple estimate gives a 
linear growth of the scattering rate with temperature. The other support for a constant (3r) with respect 
to temperature comes from the experimentally observed constant minimal conductivity, found for a wide 
range of temperatures pQ. It will be shown below that the minimal conductivity at uj — and < 00 is 
proportional to (3rj. 

Discussion of the results: I) zero temperature. With expression (jlOp . the conductivity cr 22 in Eq. (T?]) 
reads eventually 



22 ~ 4h 



-r' (? «« )+ JL)*-|( 1+ £)f 



-ui/2 



The linear increase with the scattering rate r\ is in agreement with the Drude formula 

° I + (uj/v) 2 

for large oj, except for the different power in uj. This dependence on i] is in qualitative agreement with 
the reduction of the minimal conductivity after annealing (i.e. effectively reducing 77) which was observed 
in a recent experiment by Geim and Novoselov [15] . 

The expression of the conductivity a' 22 can be studied in several limits. As a first example, in Eq. |8|) 
the limit oj — > is taken first and then the limit 77 — i> 0. Then the conductivity in Eq. Q reads 

e 2 n 2 f°° 1 1 e 2 

V 1 T kdkns~. (15) 



h TT 2 J Q (k 2 + ?7 2 ) 2 TT h 
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Figure 1: Conductivity of Dirac fermions vs. 0u> (/3 is the inverse temperature and u> the frequency) in 
units of e 2 /h. The conductivity increases with the rate f3rj = 1,2,4 (77 the scattering rate). It is assumed 
that j3rj does not depend on (3 (from Eq. (I13[) ) . 



In the next example the result in Eq. (fT 



and 



is considered. This yields for 77 10 the minimal conductivity 
S*„ Jf! f 0r (16) 
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The last result agrees reasonably well with the experimental observation of Ref. [T] if it is multiplied by 
4, the factor that is taking care of the two-fold spin and the two-fold valley degeneracy of graphene. 
II) frequency and temperature dependence. There are two asymptotic regimes with 



1 



log 



{M 2 

which implies for the conductivity 



1 + tanh 2 (/3w/4) 



1 -tanh 2 (/3w/4) 



1/8 for /3a; 
1/2/3W for f3u 





CO 



'22 



e 

8h 



0r] for (3uj ~ 

7T + 4f3r]/(3uj for /3w ~ 00 



(18) 



The result of Eq. Q14[) is reproduced when the temperature is sent to zero first. Experimentally, however, 
it is more realistic to study the DC conductivity at a nonzero temperature. Then the conductivity depends 
on the scattering rate as 0*22 x P 1 !- Remarkable is the frequency-dependent conductivity in comparison 
with the Drude formula. The Drude conductivity vanishes for large frequencies like oc w~ 2 , in contrast 



to the almost constant behavior in Eq. (fl8|) for uo < 2A and an abrupt vanishing of the conductivity if 
the frequency exceeds the energy cut-off of the Dirac fermions 2A. 

In conclusion, the Kubo formula produces a non-universal value for the minimal conductivity in 
graphene. Depending on the order of the limits, this quantity can vary over a wider range in units of 
e 2 /h. The frequency dependent conductivity of graphene at the Dirac point is exceptional, since it is 
almost constant and drops to zero when the frequency reaches the energy cut-off of the Dirac fermions. 

I am grateful to S. Mikhailov for interesting discussions and Ch. Mudry for bringing Ref. |llj to my 
attention. 
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